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INTRODUCTION 


There  are  situations  in  biochemical  research  in  which  a  standard  multi-point  equilibrium 
or  kinetic  assay  cannot  be  applied,  e.g.,  enzyme  activity  measurements  or  radioligand  binding 
studies.  This  is  particularly  true  in  cases  where  there  is  only  a  small  amount  of  tissue  available. 
To  overcome  this  problem  it  is  convenient  to  pool  regions  from  different  animals.  However,  this 
should  be  done  only  after  one  has  determined  that  there  are  no  significant  differences  among  the 
animals  in  the  study.  Any  differences  that  are  more  than  random  can  seriously  bias  the  overall 
statistical  analysis,  unless  the  investigator  is  aware  of  the  individual  variation  and  adjusts  for  it. 
As  will  be  shown,  this  statistical  analysis  involves  an  extension  of  standard  methodology.  Apart 
from  being  a  confounding  variable  in  experiments,  interest  in  understanding  the  physiological 
and  biochemical  mechanisms  that  set  apart  individual  members  of  a  species  has  been  growing 
because  analytical  tools  become  are  increasingly  sensitive  in  detecting  subtle  variance  between 
individuals. 

As  an  example  for  radioligand  binding  studies,  we  introduce  a  new  approach  using  an 
equation  that  maintains  the  simple  relationship  between  the  measured  quantities  (total  binding 
and  nonspecific  binding)  and  the  derived  variable  of  interest,  the  specific  binding  (SB).  The 
quantity  SB  is  defined  as  the  difference  between  total  binding  (TB)  and  nonspecific  binding 
(NB).  Total  binding  and  NB  are  usually  measured  with  several  replicates.  However,  there  is  no 
match  between  a  particular  TB  replicate  and  a  particular  NB  replicate.  In  fact  the  number  of 
replicates  of  the  two  quantities  may  be  different.  This  poses  no  problem  in  defining  SB  because 
only  the  mean  values  of  TB  and  NB  replicates  are  used  for  a  given  animal  i, 

[1]  SB,  =  TB"  -  NB*  (i  =  1 . n) 
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An  analysis  of  variance  (ANOVA)  could  not  be  applied  to  test  the  null-hypothesis  that 
there  are  no  statistically  significant  differences  between  SB  in  brain  sections  of  individual 
animals.  To  reemphasize,  the  original  data  do  not  have  matched  pairs  of  TB  and  NB;  therefore, 
within-group  variability  cannot  be  determined,  since  there  is  only  a  single  value  of  SB  for  each 
animal.  The  authors  suspect  that  no  simple  modification  of  ANOVA  is  appropriate  for  this 
situation.  Indeed,  here  it  is  necessary  to  formulate  this  problem  as  a  multiple  linear  regression 
model  expressed  with  the  aid  of  a  design  matrix.  This  approach,  although  possible,  is 
cumbersome  and  therefore  we  decided  to  explore  designs  beyond  the  general  linear  model. 

We  found  a  solution  for  that  particular  problem  using  a  regression  equation  with  2 
dependent  variables  (TB,  NB)  and  an  indicator  or  index  variable,  which  identifies  the  parameters 
to  be  matched  to  the  observations.  There  is  one  such  equation  for  each  animal.  Although  linear 
in  its  parameters  the  unusual  nature  of  this  regression  equation  requires  computational  flexibility 
not  found  in  standard  regression  software  packages.  However,  our  implementation  of  a  nonlinear 
regression  program  can  successfully  be  applied  to  this  and  even  more  complex  systems. 
Nonlinear  regressions  include  as  a  special  case  multiple  linear  regressions.  They  are  more 
complex  in  mathematical  theory  and  have  a  correspondingly  more  complex  computer 
implementation,  but  they  are  extremely  versatile. 

CRH,  a  41  amino-acid  peptide,  is  widely  distributed  in  the  brain  and  is  thought  to  be  the 
major  coordinator  in  behavioral,  neuroendocrine,  autonomic,  and  immunological  responses  to 
stress  (1).  Our  primary  interest  was  the  detection  of  significant  inter-individual  differences  in 
CRF  binding  in  a  relatively  small  population  of  rats  that  had  previously  been  shown  to  differ  in 
their  behavioral  response  to  centrally  administered  CRF. 
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BINDING  METHODS 


Animals.  Adult  male  Long-Evans  Hooded  rats  weighing  325-340  g  were  used 
throughout  this  study.  Animals  were  housed  in  a  light-  and  temperature-controlled  room,  with 
controlled  administration  of  food  and  water.  Stainless-steel  guide  cannulae  were  implanted  with 
access  to  the  lateral  ventricle  for  the  i.c.v.  administration  of  oCRH  (2).  The  rats  were  allowed  to 
recover  for  1  week  before  they  received  weekly  i.c.v.-injections  of  increasing  CRH  doses  (0.1  to 
10  pg)  over  a  time  period  of  2-3  months  to  determine  the  CRH  dose-response  curve  for  each 
animal. 

Materials.  Iodinated  ovine  CRH,  [125I]Tyr0-oCRH  (125I-oCRH),  with  a  specific  activity 
of  2200  Ci/mmol  was  obtained  from  New  England  Nuclear  (Boston,  MA),  and  unlabeled  ovine 
CRH  (oCRH)  and  unlabeled  rat/human  CRH  (r/hCRH)  were  purchased  from  Bachem 
(California).  All  other  standard  reagents  were  purchased  from  Sigma  (St.  Louis,  MO). 

Tissue  preparation.  Rats  were  sacrified  by  decapitation,  and  brain  regions  of  interest 
were  dissected,  weighed,  and  placed  in  20  volumes  of  ice-cold  buffer  (50  mM  Tris-HCl,  10  mM 
MgCl2,  2  mM  EGTA,  pH  7.2). 

[125I]Tyr°-oCRH  binding  assay.  The  CRH  receptor  binding  assay  was  conducted 
essentially  as  described  elsewhere  (3).  Briefly,  tissues  of  various  brain  regions  were 
homogenized  in  4  ml  of  buffer  (50  mM  Tris-HCl,  10  mM  MgCl2,  2  mM  EGTA,  pH  7.2)  using  a 
Brinkman  polytron  (setting  6  for  10  s).  The  homogenate  was  first  centrifuged  at  30,000  x  g  for 
10  min  at  4  °C;  the  resulting  pellet  was  discarded  and  this  process  was  repeated  a  second  time. 
The  pellets  were  resuspended  in  an  incubation  buffer  (50  mM  Tris-HCl,  10  mM  MgCl2,  2  mM 
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EGTA,  0.1  %  BSA,  aprotinin  (100  kallikrein  units/ml)  and  0.1  mM  bacitracin,  pH  7.2)  to  a  final 
protein  concentration  of  0.1  - 1.2  mg/ml  assay.  The  final  protein  concentration  of  each 
membrane  preparation  was  determined  using  bovine  serum  albumin  (BSA)  as  a  standard  (4). 

Incubation  assay:  100  pi  of  the  membrane  suspension,  100  pi  of  incubation  buffer 
containing  125I-oCRH  (final  concentration  0. 1  nM),  and  50  pi  of  incubation  buffer  containing 
oCRH  (5.4  nM  final  concentration)  was  added  to  1.5  ml  polypropylene  microcentrifuge  tubes 
and  incubated  in  quadruplicate  with  50  pi  of  the  incubation  buffer  or  incubated  in  triplicate  with 
50  pi  of  incubation  buffer  containing  unlabeled  r/hCRH  (1  pM  r/hCRH)  to  define  nonspecific 
binding.  All  binding  incubations  were  initiated  by  the  addition  of  membrane  protein.  The 
reaction  was  allowed  to  proceed  for  2  h  at  22  °C.  The  tissues  were  separated  from  the  incubation 
medium  by  centrifugation  in  a  microcentrifuge  for  3  min  at  12,000  x  g.  The  resulting  pellets 
were  washed  gently  with  1  ml  of  ice-cold  phosphate-buffered  saline  (PBS:  140  mM  NaCl,  2.6 
mM  KC1,  10  mM  Na2HP04,  1.76  mM  KH2P04),  pH  7.4,  containing  0.01  %  Triton  X  100.  Then 
the  contents  were  recentrifuged  for  3  min  at  12,000  x  g;  the  supematent  was  aspired  and  the 
radioactivity  of  the  pellets  were  measured  in  a  gamma  counter  at  80%  efficiency. 

Single-point  CRH-receptor  binding  was  performed  using  an  approximate  receptor- 
saturating  concentration  of  0.1  nM  [125I]-Tyr0-oCRH  plus  5.4  nM  unlabeled  oCRH  to  determine 
TB.  Nonspecific  binding  was  defined  as  the  binding  not  displaced  by  1  pM  r/hCRH.  The  total 
concentration  of  ligand  used  consisted  of  the  concentrations  of  radiolabeled  and  unlabeled  ligand, 
[L*]  and  [L],  respectively.  To  account  for  the  dilution  with  unlabeled  ligand,  SB  was  divided  by 
the  fraction  of  labeled  ligand  f:  f  =  [L*]/([L*]  +  [L]). 
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The  standard  error  of  specific  binding  (SEMsb)  represented  in  the  figures  was  calculated 
as  the  difference  between  the  means  of  two  groups  with  unequal  number  of  data  points  (TB  with 
4  replicates  and  NB  with  3  replicates).  Equation  2  is  expressed  in  terms  of  the  two  samples 
standard-deviations  (SDTb  and  SDNb)  and  the  two  sample  sizes  (njB  and  nNB) . 


[2] 


SEMsb 


=  SD 


pooled 


where  SDpooled 


/(rc TB  ~  *1 )  ‘  S D tb  +  (*\b  ~  1)-SPNB 

V  nTB  +  nNB  "  ^ 


DATA  ANALYSIS 

Our  goal  is  to  test  the  null-hypothesis  of  no  inter-individual  differences  in  specific 
binding  (SB)  of  a  given  tissue.  To  accomplish  this  we  have  formulated  two  nested  models  to 
which  we  apply  nonlinear  regression  methods.  We  chose  iterative  nonlinear  least-squares 
regression  techniques  (algorithm  according  to  Marquardt  (5)).  For  statistical  comparison  of  the 
two  models  a  partial  F-test  was  performed  that  used  certain  output-quantities  of  the  regression 
analysis.  A  modified  Scheffe-test  was  applied  to  resolve  a  significant  overall  result  (a  =  0.05, 
two-sided)  post-hoc  into  individual  components.  Statistical  power  of  the  experimental  design 
was  established  using  Monte-Carlo  randomization  techniques. 

In  general,  our  regression  method  minimizes  the  NORM,  which  is  defined  as  the 
weighted  sum  of  squared  error  (SSE)  as  in  Eq.  3. 


5 


[3] 


NORM  =  SSE(B)  =  V 

i=i  <^i 

y,  (i  =  1,  ...,n)  are  the  values  of  the  dependent  variable,  n  =  number  of  datapoints,  f  describes  a 
general  mathematical  relation  between  dependent  (y)  and  independent  (x)  variables,  B  is  the 
parameter-vector  to  be  optimized,  x*  are  the  values  of  the  independent  variable  x  and  a;  is  the 
standard-deviation  of  the  dependent  variable  at  a  datapoint  (x*,  y{).  A  non-weighted  scheme  was 
used  and  the  individual  standard-deviations  (c*)  were  set  to  1 .  An  important  note:  for  the 
regression  equation,  represents  TB;  for  some  of  the  data  and  NBi  for  other  data. 

We  selected  the  sum  of  squared  error  as  a  “goodness-of-fit”  criterion.  In  comparing 

equations  with  different  numbers  of  parameters,  the  curve  generated  by  the  more  complicated 

equation  (the  one  with  more  parameters)  will  generally  come  closer  to  the  data  points.  The 

question  is  whether  this  decrease  in  sum  of  squared  error  is  worth  the  “cost”  of  additional 

parameters  (loss  of  degrees  of  freedom).  When  comparing  hierarchical  models,  the  probability 

that  additional  parameters  are  without  effect  on  the  sum  of  squared  error  is  defined  by  an  F- 

distribution  (6).  The  F-ratio  quantifies  the  relationship  between  the  relative  increase  in  sum  of 

squared  error  and  the  relative  increase  in  degrees  of  freedom  and  is  calculated  according  to  Eq.  4. 

(SSEt  -SSE2)~ 

(dft  -df2) 

sse2 
df2 

where  SSE,  is  the  sum  of  squared  error  and  dfi  is  the  degree  of  freedom  (for  model  1  and  2,  i  = 
1,2).  The  two  values  of  SSEi  and  df;  are  obtained  from  two  analyses  to  the  same  data  set. 
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differing  only  in  the  number  of  parameters.  In  the  reduced  model  some  of  the  parameters  of  the 
full  model  are  fixed  at  zero.  If  the  reduced  model  is  correct,  the  expected  F-ratio  is  near  1.0.  If 
the  ratio  is  much  greater  than  1.0,  there  are  two  possibilities:  (1)  the  full  model  is  correct  or  (2) 
the  reduced  model  is  correct,  but  random  scatter  led  the  more  complicated  model  to  fit  better. 

The  p  value  expresses  how  rare  this  coincidence  would  be  and  answers  this  question:  if  the 
reduced  model  (Index  1)  is  really  correct,  what  is  the  chance  to  randomly  obtain  data  that  fits  the 
full  model  (Index  2)  significantly  better?  If  the  p-value  is  low,  one  concludes  that  the  full  model 
fits  the  data  significantly  better  than  the  reduced  model.  Levels  of  significance  corresponding  to 
the  calculated  value  of  F  can  be  computed  from  the  probability  function  of  the  F-distribution, 
which  has  (dfi  -  df2)  degrees  of  freedom  for  the  numerator  and  df2  for  the  denominator.  To 
compute  the  p-value  the  C-subroutine  betai  was  implemented  (7). 

In  the  present  analysis  two  non-standard  regression  models  that  are  linear  in  the 
parameters  were  defined.  One  of  them  can  be  expressed  as  a  restricted  case  of  the  more  general 
form  (nested  design).  These  models  assume  that  variances  of  TB  and  NB  are  homogeneous 
across  all  animals.  Both  models  were  fit  to  the  original  data  set,  including  TB  and  NB  points 
being  corrected  for  differences  in  protein  concentrations.  The  restricted  simple  model  consists  of 
M+l  parameters  to  be  estimated,  one  of  these  representing  NB  for  each  of  the  M  animals  and  one 
additional  parameter  for  the  common  value  SB  for  all  animals.  The  general  model  consisted  of 
2-M  parameters  to  be  estimated,  one  NB  parameter  and  one  SB  parameter  for  each  of  the  M 
animals.  To  express  the  general  model  in  a  regression  form  we  rearranged  Eqs.  2-5.  Note  that 
the  vector  B  can  represent  either  of  two  types  of  observed  binding  and  thus  Eq.  5  is  of  the 
multiple  dependent  variable  type  regression  described  above. 
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[5] 


f  NB[i]  +  k-SB 
{  NB[i]  +  k-SB[i] 


for  Model  1 
for  Model  2 


where  k 


f  0  if  datapoint  represents  nonspecific  binding 
{  1  if  datapoint  represents  total  binding 


where  the  index  i  refers  to  the  i-th  animal  (i  =  When  a  datapoint  is  a  nonspecific  binding 

measure,  k  is  set  to  zero  and  the  dependent  variable  B  represents  nonspecific  binding.  If  the 
datapoint  is  a  measure  of  total  binding  then  k  is  set  to  1  and  B  represents  total  binding.  NB  and 
SB  are  the  parameter  vectors  to  be  simultaneously  optimized. 

A  full  description  of  the  nonlinear  regression  method  may  be  found  elsewhere  (8). 

Briefly,  the  first  step  is  input  of  data  and  the  starting  values  for  the  parameters  to  be  estimated 
followed  by  calculation  of  the  model  dependent  variable  B,  as  outlined  above  (Eq.  5).  Then  the 
difference  between  the  measured  datapoints  (TB,  NB)  and  the  calculated  function  B  of  the  model 
has  to  be  squared  and  summed  over  all  data  points  (Eq.  3).  The  resulting  value  is  called  the  sum 
of  squared  error  (SSE).  In  order  to  find  the  optimal  parameters  for  the  model  function,  SSE  has 
to  become  a  minimum.  At  that  minimum,  all  partial  derivatives  from  SSE  with  respect  to  each 
parameter  becomes  zero.  Those  numerical  calculated  partial  derivatives  are  used  with  further 
operations  to  select  a  new  set  of  parameters.  Again,  the  SSE  is  found  by  taking  the  sum  over  the 
squared  difference  between  the  measured  point  (TB,  NB)  and  the  function  B.  All  this  is  repeated 
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until  relative  changes  in  the  parameter  estimates  and  successive  changes  in  SSE  are  less  than  a 
certain  tolerance  (here  10‘6),  at  which  point  convergence  and  a  final  solution  is  declared. 

The  approximate  standard  error  SEk  on  each  parameter  can  be  calculated  from  the 
diagonal  elements  of  the  covariance  matrix  COV  according  to  Eq.  6.  To  calculate  SEk,  a  derived 
quantity,  the  estimated  variance  of  the  observations  (MSE :  mean  squared  error)  has  to  be 
calculated  after  the  optimization.  SSE^  is  the  weighted  sum  of  squared  error  (Eq.  3)  at  the 
minimum  and  the  difference  between  the  number  of  measurements  and  the  number  of  parameters 
is  the  degrees  of  freedom  (df).  The  resulting  formula  is  stated  in  Eq.  6. 


[6] 


SEk  =4 


covy.  -  MSE 


SSF 

where  MSE  =  1111,1 


df 


The  parameter  coefficient  of  variation  in  percent  (%CV)  was  calculated  according  to  the 
ratio  of  the  parameter  standard-deviation  SEk  and  the  parameter  estimate  bk  times  100. 

[7]  %CVk  =  100  • 

t»k 

Correlations  in  optimization  procedures  arise  when  two  or  more  parameters  try  to  explain 
the  same  data.  The  correlation  coefficient  between  any  two  parameters  (corrki)  is  given  by  the 
scaled,  off-diagonal  elements  of  the  covariance  matrix  (Eq.  8).  The  parameters  are  mutually 
independent  when  corry  is  zero,  on  the  other  hand  if  corrid  is  ±  1,  they  are  explaining  the  data  in 
the  exact  same  way. 
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[8]  corry  = 


cov„ 


Vcov kk  -covu 


Several  statistics  calculated  during  the  optimization  procedure  may  be  used  as  indicators 
of  “goodness-of-fit.”  These  are  the  root  mean  squared  error  and  the  coefficient  of  determination. 

The  root  mean  squared  error  (RMS)  measures  the  average  deviation  around  the  fitted 
curve.  RMS  in  percent  is  defined  in  Eq.  9. 


[9]  %RMS  =  100-VMSE] 

The  coefficient  of  determination  (R2)  is  a  measure  of  the  fraction  of  the  total  variance 
accounted  for  by  the  model  and  calculated  as 


[10] 


R2  =  1- 


SSE^ 

SSYW 


where  SSE^n  is  the  weighted  sum  of  squared  error  (Eq.  3)  at  the  minimum  and  SSYW  is  defined 
as  X  [(yi  -  ymean)/  Oj]  2  with  ymean  as  the  weighted  arithmetic  mean  of  the  y-values. 

In  case  the  overall  F-test  revealed  a  p  value  lower  than  the  significance  level  (a  =  0.05, 
two-sided),  a  modified  Scheffe-test  was  performed.  Of  the  approaches  considered  for  controlling 
Type  I  error  in  multiple  comparisons,  a  number  of  tests  are  available,  but  the  Scheffe  method  is 
the  only  one  that  is  based  on  the  same  assumptions  as  the  overall  F-test  in  the  analysis  of 
variance  (c.f.  Winer  et  al.  (9)).  Regarding  our  model,  out  of  the  14  parameter  estimates  (means) 
only  7  estimates  associated  with  SB  should  be  contrasted.  The  within-group  degrees  of  freedom 
can  be  determined  as  dfwg  =  ntot  -  npar  =  49  -  14  =  35,  where  ntot  are  the  total  number  of  data 
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points  and  npar  are  the  number  of  all  parameters  to  be  estimated.  The  degrees  of  freedom  for  the 
between  group  variability  can  be  written  as  dfbg  =  M-1  =  7-1=6,  where  M  is  the  number  of 
animals.  The  F-statistic  according  to  the  ideas  of  Scheffe  can  be  calculated  using  the  SB 
parameter  estimates  and  their  associated  variances  Var(SB)  as  calculated  from  the  main  diagonal 
of  the  variance-covariance  matrix. 


[11]  F(df5s,dfws) 


(SB[i]  -  SBjjjf 
dfM(Var(SB[i]  +  Var(SB[j])) 


where  i,j  =  1,...,7  and  i  ^  j.  This  formula  assumes  parameters  SB[i]  and  SB(j]  are  uncorrelated 
(corrij  =  0)  as  is  the  case  for  our  data.  Note  that  the  variance  of  the  SB  parameter  estimates  can 
also  be  expressed  in  terms  of  the  mean  squared  error  (cf.  Eq.  6)  and  the  number  of  associated 
data  points:  Var(SB[i])  =  MSE-(l/nxB  +  1/iinb)-  The  relationship  between  the  F-distribution  and 
the  incomplete  (3-function  was  used  to  compute  a  probability  p  (7). 

In  Monte-Carlo  simulations,  random  noise  with  constant  error  was  added  to  each  point 
calculated  on  the  theoretical  mean  Btme  with  an  expected  standard-deviation  of  RMS  (c.f.  Table 
1).  The  resulting  random  datapoint  Bconst  can  be  expressed  as 


[12]  Bconst  =  B^+RMS-RND 


where 


B 


true 


f  Bnb  if  data  point  is  generated  as  NB 

[  Bnb  +  Bsb  if  data  point  is  generated  as  TB 
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Bnb  and  Bsb  are  the  optimized  parameters  of  the  full  model.  RND  is  a  normal  random  deviate 
with  a  mean  of  zero  and  a  unit  variance,  N(0,1).  Generation  of  a  uniformly  distributed  random 
number  and  conversion  into  a  normal  deviate  was  implemented  using  the  subroutines  rani  and 
gasdev  given  in  Numerical  Recipes  in  C  (7). 

All  calculations  were  performed  by  using  computer  programs  written  in  Turbo-C++ 
4.5  (Borland  Inc.)  to  run  on  a  pentium  IBM-type  Personal  Computer  under  Windows  NT. 

RESULTS  AND  DISCUSSION 

Equilibrium  single-point  radioligand  binding  experiments  were  carried  out  to  test  the 
hypothesis  of  significant  inter-individual  differences  across  animals  in  a  given  tissue  preparation. 
Six  different  brain  regions  were  chosen  and  replicates  of  TB  and  NB  were  determined.  One  data 
point  of  section  #1  was  considered  to  be  an  outlier  and  was  dropped  (more  than  3  standard- 
deviations  from  the  mean). 

A  nonlinear  regression  approach  was  chosen  to  overcome  some  implementation  problems 
under  the  general  linear  model  as  explained  in  the  section  on  Data  Analysis.  Our  approach  can 
be  formulated  and  solved  in  terms  of  two  simple  nested  regression  models  (Eq.  5).  Another 
advantage  is  the  use  of  all  original  datapoints  (TB  and  NB),  corrected  only  for  differences  in 
protein  contents. 

Starting  values  for  the  8  parameters  of  Model  1  (7  NB  and  1  SB)  were  taken  from  the 
arithmetic  means  of  the  NB  datapoints  (n  =  3)  and  from  the  arithmetic  mean  of  all  SB  across  the 
M  animals  calculated  according  to  Eq.  2  (n  =  7).  Part  of  the  starting  values  for  the  14  parameters 
(7  NB  and  7  SB)  for  the  subsequent  analysis  of  the  same  data  set  with  Model  2  were  taken  from 
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the  optimized  parameter  values  of  Model  1.  Starting  values  for  Parameters  9  to  14  (all  SB)  were 
identical  to  the  starting  value  of  Parameter  8.  All  regression  cycles  converged  after  5  iterations. 
Using  the  SSE  and  the  proper  degrees  of  freedom  of  both  models,  an  F-ratio  according  to  Eq.  4 
was  calculated  and  a  p-value  under  the  F-probability  distribution  determined.  In  case  p  <  0.05,  a 
modified  Scheffe-test  was  performed  trying  to  resolve  the  overall  result  post-hoc  for  individual 
differences. 

Certain  summary  statistics  are  presented  in  Table  1.  The  SSE  in  the  full  model  is  never 
larger  than  in  the  reduced  model.  Recall  that  RMS  is  the  square  root  taken  of  SSE  divided  by  the 
proper  degrees  of  freedom.  The  degrees  of  freedom  for  the  full  model  is  always  less  than  that  for 
the  reduced  model.  In  general  the  RMS  will  be  smaller  for  the  full  model  than  the  reduced 
model,  only  when  the  extra  parameters  significantly  reduce  the  SSE  (when  the  full  model  really 
fits  the  data  better).  As  expected,  R2  as  a  descriptive  measure  of  the  “goodness-of-fit”  is  always 
greater  for  the  full  model  than  for  the  reduced  model.  However,  the  significance  of  this 
difference  is  determined  by  the  overall  partial  F-Test  (refer  to  Table  2,  last  column).  Only  the 
two  R2  of  Section  #3  (piriformis)  reach  statistically  significant  differences  at  the  0.05  level.  The 
Section  #2  is  significant  at  the  0.1  level.  All  other  sections  show  no  significant  differences 
between  the  R2  of  both  models. 

For  each  section  the  reduced  model  has  more  precise  estimates  of  its  parameters,  reflected 
by  its  lower  %CVmax.  Even  for  the  full  model  the  estimates  were  sufficiently  precise  since 
%CV ^  was  always  less  than  2/3.  Recall,  that  CV^  is  the  maximum  value  of  all  the  ratios  of 
standard  error  of  the  parameter  divided  by  its  estimate.  As  can  be  predicted  from  a  matrix 
analysis  of  the  regression  equation,  the  standard  deviation  for  each  parameter  depends  on  the 
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sample  standard  deviation  and  number  of  observations  of  TB  and  NB  for  each  animal. 

Therefore,  in  case  all  the  animals  have  the  same  number  of  TB  (nO  and  NB  (n2)  observations,  the 
parameter  standard  deviations  for  all  TBj  were  the  same  and  likewise  for  NB;  (results  not  shown). 

In  Fig.  1  specific  CRH-receptor  densities  in  the  hypothalamus,  hippocampus,  piriformis, 
locus  coeruleus,  frontal  cortex,  and  cerebellum  of  7  individual  rats  are  compared.  Table  1 
summarizes  some  results  of  the  regression  analysis  of  the  2  models.  The  statistical  analysis  for 
model  discrimination  shown  in  Table  2  indicates  no  significant  differences  between  individuals 
at  the  0.01  level.  However,  at  the  0.05  level  the  piriformis  region  reveals  significant  inter¬ 
individual  differences,  and  the  hippocampus  region  is  significant  at  the  0.1  level.  For  the 
piriformis  section,  a  further  post-hoc  contrast  among  all  rats  was  performed  using  a  modified 
version  of  the  Scheffe-test  (c.f.  Section  Data  Analysis).  In  Table  3  the  results  for  sections 
piriformis  are  presented. 

Even  if  important  differences  are  there,  one  might  not  obtain  statistically  significant 
differences  in  the  experiment.  Just  by  chance,  the  data  may  yield  a  p- value  greater  than  the 
significance  level  a.  When  interpreting  a  result  of  an  experiment  that  found  no  significant 
difference,  one  must  know  how  much  power  the  study  had  to  find  various  hypothetical 
differences,  if  they  existed.  The  power  depends  mainly  on  the  sample  size  and  the  amount  of 
variation  within  the  groups,  as  quantified  by  the  standard-deviation.  One  way  to  determine 
power  is  to  analyze  many  experiments  that  have  the  original  number  of  sample  size  (nt0t),  but 
different  artificial  values,  and  calculate  a  p  value  for  each  experiment.  Statistical  power  is  the 
fraction  of  these  experiments  that  have  a  p  value  less  than  a  and  thus  are  declared  statistically 
significant.  Table  2  shows  the  result  of  empirically  determined  power  for  the  0.05  significance 
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levels.  It  is  evident  that  only  the  data  from  Section  3  (piriform  cortex)  have  sufficient  power  to 
detect  a  possible  true  difference  in  an  analysis.  However,  contrasting  the  individual  means  of 
this  brain  section  do  not  indicate  which  is  different  (Table  3). 

In  summary,  statistical  analysis  of  possible  inter-individual  differences  in  SB  between  rats 
failed  to  reject  the  Null-hypothesis  at  the  0.01  level.  However,  in  one  region  (piriformis  cortex) 
the  Null-hypothesis  could  be  rejected  at  the  0.05  level  and  in  another  region  (hippocampus)  at  the 
0.1  level  (see  also  Table  3).  A  further  multiple  comparison  of  the  section  piriformis  could  not  be 
resolved  for  individuals.  Results  of  a  statistical  power-analysis  showed  that  a  high  p-value  for 
the  overall  partial  F-test  is  accompanied  by  a  fairly  low  power.  This  indicates  that  if  in  fact  a 
difference  would  exist,  the  chosen  experimental  design  would  not  be  able  to  detect  it  most  of  the 
time.  In  such  instances,  an  application  of  this  new  approach  may  be  of  use  in  providing 
appropriate  analytical  methodology.  The  application  of  techniques  that  lend  themselves  to  the 
analysis  of  individual  differences  in  animals  will  increase  appreciation  for  how  these  differences 
may  be  elucidated  at  a  molecular  and  biochemical  level  of  an  organism. 
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Table  1 

Comparison  of  Regression  Results  between  the  Two  Models 
throughout  Different  Brain  Sections. 

# 

Section 

Model  1 

Model  2 

%CVmax 

%RMS 

R2 

%cvmax 

%RMS 

R2 

1 

hypo 

9 

6.1 

0.81 

42 

5.8 

0.94 

2 

hippo 

12 

8.7 

0.92 

62 

8.0 

0.94 

3 

piri 

7 

8.7 

0.93 

16 

7.7 

0.96 

4 

lc 

10 

8.9 

0.81 

37 

9.4 

0.82 

5 

fctx 

5 

6.7 

0.95 

11 

6.3 

0.96 

6 

cere 

9 

13.3 

0.83 

34 

12.9 

0.86 

Model  1  =  restricted  Model  (8  parameter ) ;  Model  2  =  full  model  (14  parameters );  °/ = 
maximum  value  coefficient  of  variation  as  calculated  according  to  Eq.  7;  %RMS  =  root  mean 
square  in  percent  as  calculated  using  Eq.  9;  R2  =  coefficient  of  determination  as  computed  using 
Eq.  10 

Abbreviations:  1.  hypo  =  hypothalamus;  2.  hippo  =  hippocampus;  3.  piri  =  piriformis;  4.  Ic  = 
locus  coeruleus;  5.  fctx  =  frontal  cortex;  6.  cere  =  cerebellum. 

For  brain  sections  #  1-6  the  maximum  correlation  coefficient  between  the  parameter  estimates  of 
model  1  was  -0.40  and  -0. 76  for  Model  2. 
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Table  2.  Results  of  Significance  Tests  for  Inter-Individual  Variability 
in  Different  Brain  Sections  of  Rats. 

# 

Section 

SSE 1 

dfi 

F(dfrdf2,df2) 

P 

Power 

sse2 

df2 

a=  0.05 

1 

Hypo 

1.394 

40 

1.758 

0.138 

0.56 

1.064 

34 

H 

Hippo 

2.801 

41 

2.255 

0.061 

0.70 

11 

2.020 

35 

B 

Piri 

5.473 

41 

2.823 

0.024* 

0.82 

■ 

3.688 

35 

H 

Lc 

2.073 

41 

0.246 

0.958 

0.11 

■ 

1.789 

35 

A 

Fctx 

3.703 

41 

1.947 

0.100 

0.64 

1  S 

2.777 

35 

Q 

Cere 

10.782 

41 

1.441 

0.227 

0.46 

m 

8.647 

35 

SSEj  =  sum  of  squared  error  of  Model  1  ( reduced  model)  at  the  minimum;  SSE2  =  sum  of 
squared  error  of  Model  2  (full  model)  at  the  minimum;  dfi  and  df2  are  the  respective  degrees  of 
freedom.  The  F -ratio  is  revealed  according  to  Eq.  4.  Statistical  power  expressed  as  fraction  is 
empirically  determined  with  2000  independent  Monte-Carlo  simulations  (a  two-sided). 
Abbreviations:  1.  hypo  =  hypothalamus,  2.  hippo  =  hippocampus,  3.  piri  =  piriformis,  4.  Ic  = 
locus  coeruleus,  5.  fctx  =  frontal  cortex,  6.  cere  =  cerebellum. 

*p  <0.05 
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T 

able  3.  Results  of  a  modified  Scheffe-Test  for  Inter-Individual 

Variability 

in  Area  Piriformis  in  Brains  of  Different  Rats. 

RAT # 

1 

2 

3 

5 

6 

7 

MEAN 

122.9 

134.6 

130.5 

117.2 

75.9 

98.8 

1 

122.9 

1 

1.00 

1.00 

2 

100.9 

0.95 

1 

EH 

0.82 

0.99 

0.91 

1.00 

3 

134.6 

0.99 

1 

1.00 

0.98 

0.11 

0.65 

4 

130.5 

1.00 

0.82 

1.00 

1 

1.00 

0.17 

0.77 

5 

117.2 

1.00 

0.99 

0.98 

1.00 

1 

0.49 

0.98 

6 

75.9 

0.33 

0.91 

0.11 

0.17 

0.49 

1 

0.94 

7 

98.8 

0.92 

1.00 

0.65 

0.77 

0.98 

0.94 

1 

Multiple  comparison  results  using  a  modified  Scheff e-test  comparing  animals  after  the  overall  F- 
test  revealed  a  probability  p  <0.05  (cf  Table  2).  The  F -ratio  was  calculated  according  to  Eq.  4. 
The  shaded  area  provides  information  about  a  certain  rat  with  its  estimated  specific  binding 
( mean  value),  expressed  infmol/mg  protein.  The  white  area  shows  the  computed  probability  of 
testing  individual  differences  (c.f  Eq.  11). 
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Figure  1.  Individual  levels  of  specific  CRH-receptor  sites  for  [125I]-Tyr0  oCRH  in  rat  brain 
hypothalamus,  hippocampus,  piriformis,  locus  coeruleus,  frontal  cortex  and  cerebellum  after 
repeated  treatment  with  i.c.v.  oCRH,  expressed  as  mean  and  SEM.  Membrane-suspensions  of 
each  section  were  incubated  for  2  h  at  room  temperature  in  the  presence  of  0.1  nM  [125I]-Tyr0 
oCRH  and  5.4  nM  unlabeled  oCRH  (mean  protein  content  ranged  from  0.1  to  0.3  mg  protein/ml 
assay).  Abscissa:  seven  individual  rats.  Ordinate:  specific  binding  in  fmol/mg  protein  as  revealed 
from  equation  1.  SEM  was  calculated  according  to  Eq.  2. 
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